Negative refraction in nonlinear wave systems 
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Abstract 

People have been familiar with the phenomenon of wave refraction for several centuries. Recently, 
a novel type of refraction, i.e., negative refraction, where both incident and refractory lines locate 
on the same side of the normal line, has been predicted and realized in the context of linear optics 
in the presence of both right- and left-handed materials. In this work, we reveal, by theoretical 
prediction and numerical verification, negative refraction in nonlinear oscillatory systems. We 
demonstrate that unlike what happens in linear optics, negative refraction of nonlinear waves does 
not depend on the presence of the special left-handed material, but depends on suitable physical 
condition. Namely, this phenomenon can be observed in wide range of oscillatory media under 
the Hopf bifurcation condition. The complex Ginzburg-Landau equation and a chemical reaction- 
diffusion model are used to demonstrate the feasibility of this nonlinear negative refraction behavior 
in practice. 
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It is well known that light can refract at interface of two media. For several centuries only 
positive refraction has been found in nature, where incident line and refractory line locate on 
either side of the normal line of the interface. Very recently, a new type of optical refraction- 
negative refraction, where both incident and refractory lines locate on the same side of the 
normal line of the interface, has been theoretically predicted [l] and experimentally observed 
[2I when refraction occurs between a right-handed material and a left-handed material where 
optical waves propagate inward to the wave source. In addition, some new observations have 
been reported in nonlinear wave propagation, and one particularly interesting phenomenon 
is antispiral waves |3|, In contrast to spirals (which are familiar to people) 

where waves propagate outwards from the tips, antispiral waves propagate inward to the 
spiral tip, the wave source. Because the existence of inwardly propagating waves is key 
to generate negative refraction in linear optics, we expect that the novel phenomenon of 
negative refraction may also be observed in other nonlinear wave systems where antispirals 
and inwardly propagating waves can be found. 

In this letter, we go further from the above phenomena of negative refraction in linear 
optics and antispirals in nonlinear waves propagation by revealing a new type of negative re- 
fraction in nonlinear oscillatory systems. We demonstrate that unlike what happens in linear 
optics, negative refraction of nonlinear waves does not require special left-handed materials, 
rather, it can be found in wide range of media under suitable nonlinear condition of Hopf 
bifurcation. To demonstrate this new phenomenon, we use the complex Ginzburg-Landau 
equation (CGLE) as our example. With this model, we are able to freely generate, by pe- 
riodic forcing, both normal waves (NWs) and antiwaves (AWs), of which NWs propagate 
outward from the wave source while AWs propagate towards the wave source. An interesting 
feature of negative refraction is observed when waves on either side of the interface have 
different propagation characteristics, one AW and the other NW. Moreover, we use a typ- 
ical chemical reaction diffusion model, the Brusselator, to verify the prediction of negative 
refraction of CGLE. This suggests possible experimental realization of negative refraction 
in nonlinear chemical reaction diffusion systems. 



We consider 2D CGLE with local periodic forcing 
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where no-flux (Neumann) condition is applied to all boundaries. The last term on the right 
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side of the equation, an additive periodic forcing with constant amplitude F, serves as a 
wave source applied to the x = stripe of the medium. F needs to be sufficiently large 
so that periodic waves can be stimulated, which will then propagate as planar waves into 
the unforced region of x > 0. We use F — 1 throughout this paper. The forcing frequency 
uji n has its sign, which denotes the rotating direction of the forcing in the complex plane, 
"+" clockwise and "-" anticlockwise. In the absence of the source term [F = 0), Eq.(l) 
describes the onset of homogeneous oscillation through Hopf bifurcation, and this equation 
is used extensively as a typical model for descri 
propagating waves and pattern formations flO . 

In Fig.l, we fix (3 = —1.4 and demonstrate distinctive dynamical behaviors of the systems 
for different a and uj in . Some parameter combinations generate NW (Fig. 1(a)), while others 
generate AW (Fig. 1(b)). The parameter regions in Fig. 1(c) indicated by NW and AW 
show NW and AW outputs, respectively, with co = u in (u denotes the output frequency). In 
the shadowed region, the periodic forcing cannot control the system, and thus the system 
motion under periodic forcing is either periodic with oj ^ uji n , or chaotic with broad band 
frequency spectrum. In this study we are not interested in the shadowed region, and will 
focus on the NW and AW regions only where u = u in . Interestingly, all the NW and AW 
regions locate in the vicinity of the diagnol line u = Ui n = a. Since a is the bulk oscillation 
frequency (u = a) of the local dynamics, all the NW and AW regions of Fig. 1(c) satisfy 
resonance condition (iOi n ~ u = a). Moreover, in Fig. 1(c) it is clear that AWs appear only 
under the condition 

uujo > and |a>| < \ujq\ = \a\ (2) 

NWs can be observed under the conditions of ujujq < or |o;| > \ujq\ for ujujq > 0. The 
boundaries between the parameter domains of NW and AW can be determined by the sign 
change of the phase velocity [z| v p h aS e — f ■ The transitions between NW and AW thus 
occur at two boundaries: u — and k — 0. The former boundary justifies the transition 
between AW and NW with stationary Turing patterns uj = u in = at the transition, and 
the latter defines the transition between AW and NW with homogeneous oscillation (k = 0) 
at the transition, which shrinks to a single point of a — /3 — Ui n in Fig. 1(c). In summary, 
we find that, with local forcing method, we can generate NWs and AWs freely. For instance, 
given an oscillatory medium with fixed parameters, we can produce both NW and AW by 
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FIG. 1: In this and all the following figures, we apply no-flux (Neumann) boundary condition. 
(3 = —1.4, F = 1.0 are used in Figs.l and 3 and 256 x 256 rectangular geometry is used for 
numerical simulations in Fig.l. (a & b) Numerical results of spatiotemporal patterns at y = 128, 
x G [0, 256] with different periodic forces applied on the left boundary of the rectangular system. 
Arrows denote the wave propagation directions. Space descritization of 256 x 256 pixes and time 
step At = 0.05 are used for numerical simulations, (a): Normal forward propagating waves (NWs). 
a = 0.2, uj = Ui n = —0.15. (b): Planar antiwaves whose phase propagates toward the wave source. 
a = 0.2, u = ijJi n = 0.15. (c): Distributions of different characteristic regions in Ui n — a plane. 
Shadowed region: no stable waves with to = uji n can be realized in Eq.(l) whatever the input 
frequency U{ n . NW white region: normal waves of u = Wj n can be observed; AW white region: 
antiwaves of u = u>i n can be observed. 

changing the forcing frequency u>i n only. This is different from antispirals in autonomous 
systems which produce either NW or AW (not both) with a single spiral frequency selected 
at the given set of parameters. 

Now we come to the problem of nonlinear wave refraction. We consider planar waves 
passing an interface of two media with different parameters a, f3 (Fig. 2). Analogous to 




negative refraction 



positive refraction 



FIG. 2: Schematic figure of nonlinear refraction. The left and right regions are CGLE systems 
with parameters a%, Pi and fo, respectively. The middle line from left below to right up is the 
interface. The line marked with k\,k2, (fo) are incident line, positive refractory line, and negative 
refractory line, respectively. 

optics, the incident line and refractory line are perpendicular to the wave fronts in the two 
media. The incident angle (9i) is denned as the angle between the normal line and the 
incident line, which is positive definite. The refractory angle (62) is defined as the angle 
between the normal line and the refractory line, which is positive if the refractory and 
incident lines locate on either side of the normal line, and is negative otherwise. 

In linear optics, the light refraction relation can be precisely computed, based on the 
continuity conditions at the interface. For nonlinear systems, there is no general theory to 
describe the law of nonlinear wave refraction although the continuity conditions are still 
valid at the interface. The difficulty is that no analytical solutions can be obtained around 
the interface of two nonlinear media. We approximate the solutions of Eq.(l) as planar wave 
solutions 1 161 



in the left and right regions of Fig. 2, respectively. The actual solutions around the interface 
certainly deviate from these planar wave solutions. For instance, the amplitude of \A\\ and 
\A 2 \ must change considerably in the interface region for satisfying the continuity conditions. 
Our assumption is that the interface causes only very weak influence to the phase distribu- 
tions of the planar waves (this assumption is fully confirmed by our numerical simulations). 




(3a) 
(3b) 
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This assumption allows us to predict the nonlinear wave refraction law by applying the phase 
continuity hypothesis of the planar wave solutions on the interface: 

— u)\t + k\ ■ f — -u 2 t + k 2 -f (4) 

Here, we consider the case of a>i = u 2 only, otherwise the refractory pattern cannot be 
stationary \v\ . Therefore, we have k\ - r = k 2 - r which leads to k\ sin#i = A^sin^- 

There are several interesting conclusions we can draw from Eq.(4). First, with the sign 
of uj fixed (e.g., positive u>) k x and k 2 can be positive (NW) or negative (AW). It tells us 
that if ki and k 2 have the same sign positive refraction will be observed. If k\ and k 2 
have different signs, i.e., from NW at left to AW at right or vise versa, nonlinear negative 
refraction, an analog of propagation of light from right- to left-handed material (or vice 
versa), emerges. Second, with the help of dispersion relation, we can define from Eq.(3b) an 
effective refractory index similar to that in optics as 

«H = ±n\/| EZ («) 

\u)\ V p — a 

where "+" for NW (ww < or \u\ > \u \ = \a\ and uju q > 0) and "-" for AW (ujuj > and 
| a; | < \ujq\). The generalized refraction law for nonlinear waves in CGLE can be expressed 
as 

sin 6>i _ ria ^ 
sin 9 2 n\ 

In the following, we conduct numerical simulations of (1) to test all the above analytical 
results. In our simulations, we divide a rectangular region into two parts like Fig. 2. In 
Fig. 3(a), a pattern of positive refraction of nonlinear waves is illustrated clearly with direct 
numerical simulation of Eq.(l). In Fig. 3(b), we plot the relation of incident and refractory 
angles in Q\ — 8 2 plane. The circles are numerical data and the solid line is theoretical 
prediction, which perfectly coincide with each other, fully confirming the assumption of 
Eq.(4) and the prediction of (6). 

Now we demonstrate negative refraction problem by simulation. In Fig. 3(c) and 3(d), 
we take parameter combinations at which the incident waves are NWs while the refractory 
waves are AWs. Fig. 3(c) shows the pattern of negative refraction analogous to the novel 
negative refraction in linear optics [2j. Unlike the optical case, here the negative refraction 
is observed with nonlinear waves. In Fig. 3(d), we perform the same computation as Fig. 
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FIG. 3: (a & b) Positive refraction, (a): Wj n = 0.2, 0\ = 0.1. Pattern of positive refraction with 
a\ = 0.5, j3\ = 0.0 (AW medium) and 02 = 0.5, fa = —1-8 (AW medium). The red line is the 
interface of two media. The blue line is the normal line perpendicular to the interface. The other 
two lines marked with arrows are the incident and the refractory lines, respectively. The arrows 
denote the directions of phase velocities, (b): The relation between the incident angle (#1) and the 
refractory angle (62)- a\ = 0.5, 0\ = —1.0, «2 = 0.5, = —1-4, and u>i n = —0.02. Circles represent 
numerical results and the solid line shows theoretical predictions of Eq.(6). Both results coincide 
with each other perfectly. (c)&(d) u)i n = 0.02. Negative refraction with a\ = —0.5, j3\ = 1.4 (NW 
medium) and oli = 0.5, P2 = —1-4 (AW medium), (c): Pattern of negative refraction with 0\ = 0.5. 
(d): The same as (b) with negative refraction considered. 

3(b) by considering negative refraction, and numerical results are in perfect agreement with 
the theoretical predictions for negative refraction. 

Since CGLE appears universally in nonlinear extended systems around Hopf bifurca- 
tion from homogenous stationary states to homogeneous oscillations, the negative refraction 
phenomenon found in CGLE is expected to be observed experimentally in wide range of 
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practical nonlinear and nonequilibrium systems. Chemical reaction diffusion (RD) systems 
are the most probable candidates of this kind. Here we take a well known chemical reaction 
model, the 2D Brusselator, as our example Q 



In these equations, a and b are two constant chemical components, and u and v chemical 
variable components. All a, b, u and v are positive definite (physical constraint). Without 
periodic forcing, this system undergoes Hopf bifurcation to homogeneous oscillation at the 
parameters b = 1 + a 2 . At the Hopf instability, the frequency of the oscillation reads ujh = a. 
We simulate Eq.(7) with a parameter set slightly beyond Hopf bifurcation, and find positive 
refraction phenomenon shown in Fig. 4(a). In Fig. 4(b), we use ur = 1.05 and parameter 
sets of a = 1.0, D = 0.5 and a = 1.2, D = 2.25 for the two chemical media on either side 
of the interface, respectively. The chemical parameter set can be transformed to CGLE 
parameter set by a = a-§^j-, (3 = 7 f a ^2+a^) anc ^ Uin = £ = b - (I + a 2 ) [18j. From this 

correspondence, we justify NW and AW for the left and right media of Fig. 4(b), respectively, 
and expect to observe negative refraction. Indeed, we do reveal this interesting phenomenon 
in the simulation of this nonlinear RD system in Fig. 4(b). 

In conclusion, we have generated antiwaves in CGLE and reaction diffusion systems by 
local periodic forcing. At fixed parameters, we can freely produce normal or anti waves by 
adjusting the input frequency uj in . A frequency dependent refraction index n(u) is defined 
for CGLE, which takes positive and negative values for NW and AW, respectively. When 
waves propagate from positive to negative n(u) or vice versa, novel negative refraction 
emerges. The negative refraction phenomenon predicted in linear optics depends on the 
existence of left-handed materials which have not yet found so far in nature. In our case, 
nonlinear negative refraction can be observed in wide range of natural oscillatory media 
under suitable physical condition, i.e., Hopf bifurcation and resonance conditions. Since the 
complex Ginzburg-Landau equation and nonlinear waves are common in nature, we predict 
that the negative refraction phenomenon explored in this letter can be widely observed, and 
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FIG. 4: Numerical simulation of Eq. (7), the periodically forced Brusselator. Time step At = 0.01. 
400 x 300 physical space is descritized with 400 x 300 pixes and the forcing is applied on the sites 
of the left line of the rectangular space, b = 3.2. (a) Positive refraction observed between two AW 
media, a = 1.2, D = 1.7 for the left medium and a = 1.2, D = 2.5 for the right one. lor = 1.07, 
9\ = 0.2. (b) Negative refraction observed between NW (left) and AW (right) media, a = 1.0, 
D = 0.5 for the left medium and a = 1.2, D = 2.25 for the right one. ujr = 1.05, 9\ = 0.3. 

that it can be realized in experiments, such as in chemical reaction experiments. 
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